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Простая аналитическая формула решения 
уравнений для двухмерного 
разностно-дальномерного метода 
определения координат 


В противоположность широко известным итерационным методам решения уравнений для двухмерного 
разностно-дальномерного метода определения координат предлагается простая аналитическая формула 
решения. 


Введение 


Разностно-дальномерный метод определения координат широко применяется в 
навигации и системах мониторинга источников радиоизлучения. Он имеет однознач- 
ное решение при использовании данных о времени прихода сигнала при четырех (или 
более) позициях измерения. В настоящее время известны многочисленные итерацион- 
ные способы решения уравнений этого метода [1-3]. 

Между тем эти уравнения имеют простое аналитическое решение, которым ав- 
тор пользуется на протяжении ряда лет. 

Целью данной работы является получение простой аналитической формулы для 
вычислений в устройствах, требующих быстродействия и имеющих ограниченный вы- 
числительный ресурс. 


Вывод формулы решения 


Особенно просто это решение выглядит в случае двумерной задачи. Такая задача 
до сих пор актуальна в автономных системах ближней морской навигации, в системах 
наземного мониторинга источников радиоизлучения, в радиоразведке и в других прак- 
тических задачах. 


Пусть { жж хьУЬх,, у, },{х, ‚ У, } есть априори известные координаты со- 
ответственно нулевой, первой, второй, третьей позиции измерения, а {х,, у; } — коор- 


динаты цели. Координаты целей подлежат определению по измеренным величинам 
задержек 11, 12, 13 времени между сигналами, которые прошли по пути по ломаной 
{ху ху} {ху} а = 1,2,3), и сигналами, прошедшими по прямому пути 
{ху} - {х › У #. В некоторых случаях измеряются разности #, прихода сигналов меж- 
ду пунктами {х,, у} и{х,, у, }. Задачи эквивалентны между собой, так как между Е 
и г, имеется связь: 


РОТ НИ: (1) 
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где с — скорость света ( в дальнейшем принимаем с=1), 
т, расстояние между пунктами { ху фи{х, у} Традиционные способы реше- 


ния задачи в конечном счете сводятся к следующему. Если обозначить У={х,, у, } 
вектор положения цели, а А, (1=1,2,3) матрицу квадратичной формы у * 4, *У=1, по- 
рождающую гиперболу 


о 
А = 1 |? (2) 
"в 
где а, — большая полуось гиперболы, равная 
ГА 
а, = 5 (3) 
Ь, — малая полуось гиперболы, такая, что 
Ы=с-а,, (4) 
где 
и 
бы 5 , (5) 


то решение навигационной задачи обычно сводится к отысканию общего решения сис- 
темы трех уравнений из квадратичных форм 


(у -7/2) *Т*А*Т*(-7/2)=1 (=1,2,3), (6) 
где 
| Хх; .| (1) 
У Хх, 


Хх. 
=" | (8) 
У; 


Вообще, без уменьшения общности всегда можно положить { х,, у, }={0,0}. 


Уравнения (6) представляют собой систему трех квадратичных уравнений. Ее реше- 
ние различными итерационными методами не составляет больших трудностей, однако 
значительное количество циклов итераций может привести к определенным затратам 
времени. 

Получим формулу, по которой вместо итерационного решения (6) можно полу- 
чить прямое решение непосредственно. 

Обозначая 


ВЕ Ут}, В =А-Я, п", =./(,5,), В, = вв) ‚В = (В, В) (9) 


(здесь скобки подразумевают скалярные произведения) из уравнения разности хода сиг- 


налов 
ве (10) 
С учетом (1) при с=1 имеем 
ВЮ: (0 
Возводя (11) в квадрат, получаем 
В? = +7 +2*А*Е. (12) 
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С другой стороны, из уравнения треугольника {0,0} — {х‚,у, } — {х,, у, } имеем 


ЖЕ —2*1**(Е-6,), 
где 
= (Ру, * Лт, 
ег = (х; У, *1)/Е 
есть единичные вектора, такие, что 
(е,,е,) =1, (2.,е,) =1. 


Введем обозначения 


Тогда 
(2) =(ж*х+у,*У)/н. 
Вычитая из (12) уравнение (13), получаем три уравнения 
2 2 
= „НЫ 
2 (ичк*(еье,)) 


(=1.2,3). 


Вводя обозначение 


Е? 2 
иг Е 


1-3 


(13) 


(14) 
(15) 


(16) 


(17) 


(18) 


(19) 


(20) 


и попарно вычитая величины обратные А из уравнений (19), получим три уравнения 


(ху, *У)/ и, -(х,*х+у, *Уу и, = Ви, -Е/и, (УЕ, 3 1 У), 


которые лучше выглядят как 
Ах *у=С,, (12,312), 
где 
А; =и,*х,-и,*х,, 


уж Щи 
В,, =и, Уи тУ,, 


1, 


= * > * 
брате. 


(21) 


(22) 


(23) 
(24) 
(25) 


Нетрудно показать из (6), что ранг системы (22) равен 2 (причем даже в трехмер- 
ном случае). Таким образом, чтобы найти х и у, достаточно в (22) решить любую па- 


ру уравнений. Решение (22) находим по правилу Крамера. 


= Аха 
и 
ь: ДУ 
А 
где 
А. = А, В, 
713 › 
, А, В | 
С, Г В 7 
Ах и: у у 5 
С В й 
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(28) 
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АС 
Ду = ры (30) 
Здесь везде каждое из 1, /,К может принимать значения 1,2,3, причем 1= ], 1*К, 
1=К. 
Следующим шагом находим А, подставляя (26) и (27) в (18) и (19). 
И окончательно, учитывая, что А = №* ё,, находим х, и у; 
х. =В*х, (31) 
уг =В*У. (32) 


Проведя подстановку (26) и (27) в (31) и (32), и вводя обозначения (здесь | | об- 


означают детерминант) 


Ех 
Е Хх. 
т, = Ее: (33) 
У 
г, я 
Е и, 
1; и; 
п, = вЫ (34) 
21 У, 
г, У; 
после несложных преобразований, получаем окончательно 
п.п. 
х= а. | (35) 
т, — Тя 
У, =и,-т, *х,. (36) 


Итак, алгоритм решения таков. 
Шаг 1. Зная текущие позиции измерителей и измерив задержки т,, (1=1,2,3), вычисля- 
ем Е ЕТ. - к и затем по (20) и, =2*и*т,-т2. 

В случае, когда #, известно из измерений, по (1) и (20) вычисляем и, =2***г, —г. 


Шаг 2. Вычисляем любую пару отношений детерминантов по (32) и по (33). 
Можно выбрать лучшую из соображений точности вычислений. 

Шаг 3. Вычислить искомые координаты цели по (34) и (35) 
Например: 


х,=—в в (37) 


(38) 
В заключение отметим, что уравнения (20) справедливы и в трехмерном случае. 
Здесь так же достаточно найти только два неизвестных х и у, так как третья ко- 


= — * 
Уп бт ох. 


ордината 2 единичного вектора связана с ними уравнением х’ + у’ +2’ =1. Однако 


вид решения получается довольно громоздким. 


64 «Искусственный интеллект» 172011 


Простая аналитическая формула решения уравнений... 1-3 


Выводы 


Предложенная формула может быть использована для вычислений в устройствах, 
требующих быстродействия и имеющих ограниченный вычислительный ресурс, а так- 
же как хорошее начальное приближение для итерационных методов. 
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В.М. Зуев 

Проста аналтична формула вир!шення р!внянь для двомфного 

р1зницево-далеком!рного методу визначення координат 

У протилежнеть широко в1домим 1теращйним методам виритення равнянь для двомпрного р1зницево- 
далекомрного методу визначення координат пропонуеться проста аналтична формула рипення. 
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